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Abstract 

A number of important results of studying large deformations of hyper-elastic shells 
are obtained using discrete methods of mathematical physics In the present 

paper, using the variational method for solving nonlinear boundary problems of statics 
of hyper-elastic membranes under the regular hydrostatic load, we investigate pecu- 
liarities of deformation of a circular membrane whose mechanical characteristics are 
described by the Bidermann-type elastic potential. We develop an algorithm for solv- 
ing a singular perturbation of nonlinear problem for the case of membrane loaded 
by heavy liquid. This algorithm enables us to obtain approximate solutions both in 
the presence of boundary layer and without it. The class of admissible functions, on 
which the variational method is realized, is chosen with account of the structure of for- 
mal asymptotic expansion of solutions of the corresponding linearized equations that 
have singularities in a small parameter at higher derivatives and in the independent 
variable. We give examples of calculations that illustrate possibilities of the method 
suggested for solving the problem under consideration. 

1 Definition of the strained deformed state of a circular 
membrane under the regular hydrostatic load 

Let us consider a circular membrane of the radius i?o of incompressible, isotropic and 
hyperelastic material having the small constant width Hq. Let a hydrostatic load Q be 
applied to the membrane with an inflexibly fixed contour. To describe the geometry of 
the deformed membrane, let us introduce the cylindrical coordinate system Ozrjr, the axis 
Oz of which coincides with the symmetry axis of the membrane. It follows from the axial 
symmetry of the problem that principal deformation directions at any point will coincide 
with meridians, parallels and normals to the deformed surface and all the parameters 
of the strained deformed state will be functions of initial distances s of points of the 
membrane from its symmetry axis only. Denote the principal degrees of lengthenings in 
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these directions by Ai, A2 and A3, respectively. Then, according to the nonhnear theory 
of elastic membrane Q, interior stresses Ti and T2 of the deformed shell in the direction 
of its meridian and parallel can be determined from the relations 

T. = 2/.0A3 (A? - Ai) + AL,|^) i^ = 1, 2), 

, r{s) , 1 h{s) (1.1) 




S A1A2 Hq 

/i=A? + Ai + Ai, h = + X^^ + X^^ W = W{h,l2). 

Here z{s) and r{s) are axial and radial components of the radius of the generatix vector 
of the deformed shell, is a deformation energy function for the material of the shell, Ii 
and I2 are deformation invariants, and h{s) is the width of the membrane in the deformed 
state. 

Let us choose the deformation energy function W{Ii,l2) in the form of its four-term 
approximation proposed by V.L. Biderman Q, i.e., 

W = Ci{h - 3) + C2{l2 - 3) + Giih - 3)' + C^ih - 3)3, (1.2) 

where Ci (i = 1, . . . ,4) are physical constants that are determined experimentally. 

By using function (|1.2| ) one can sufficiently precisely describe the law of deforming for 
a number of types of elastic materials up to their break. 

If the origin of the coordinate system coincides with a pole of the deformed membrane, 
then the hydrostatic pressure acting on the shell, is of the form 

Q = C* - Dz, C* = P* - P°, D = pg. (1.3) 

Here p is the density of liquid, g is the acceleration of gravity, P* is a constant component 
of pressure in liquid, and is the pressure of gas onto the shell. 
Let us introduce the following dimensionless parameters: 

K,r*,z*} = {s,r,z}/i?o, T* =T,/{2Ciha), 
W* = W/Cu Q* = QRo/{2Ciho). 

In what follows, we shall use these parameters omitting asterisks for the sake of sim- 
plicity. 

Movements corresponding to the static equilibrium state are singled out from all ad- 
missible movements by the fact that stationary values of the functional [§] 



Wih,l2)s-Qr'^ 
as 



ds (1.4) 



correspond to them. 

Having placed the origin of the coordinate system Ozr]r at the center of the undeformed 
membrane, we present the hydrostatic load as 



Q = C-Dz, 
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where difference between the constants C and C* in formula (|l.3| ) is due to the shift of 
the origin of the coordinate system. In this case, we should seek a stationary value of 
functional (|1.4|) within the set of functions satisfying the boundary conditions 



dz 

ds 



0, 



z{l)=0, r(0)=0, r(l) = 1. 



(1.5) 



s=0 



Utilizing usual tools of the calculus of variations, one can show that the Euler equations 
for functional (|1.4|) are: 



(1.6) 



(i dv 

— (rTi) = T2 — , hTi + k2T2 = Q. 
ds ds 



Here ki and /c2 are principle curvatures of the deformed median surface of the membrane 
calculated by the formulae 



/ d'^r dz dr d^z 



- 



\ ds"^ ds ds ds"^ 



A 



1 ' 



dz 
ds 



Following Rietz to find an extremum of functional (1.4) we present functions z[s) and 
r(s) in the form of the following expansions: 



'{s) = ^ XkUk{s), r{s) = S + ^ Xk+mVk{s), 



(1.7) 

k=l k=l 

where in view of ( |1.5D , the coordinate functions Ufc(s) and Vk{s) have to obey the conditions 

n^O) = Uk{l) = Vk{Q) = Vk{l) = 0. (1.8) 

Then constants Xk (A; = 1, 2, ... , 2m) forming 2m-component vector x will be determined 
from the nonlinear algebraic system 



g{x) = 0. 

In this case, components of the vector-function g are of the form 



(1.9) 



9i+m 



dz dui dr 

U{Xi, \2)-: — 1 QA2— 

ds ds ds 



dr dvk 



s ds, 



ds ds 



A, 



dz 



^^(Ai, A2) — + C/(A2, Ai)^ + gA2— V. 



ds 



sds. 



(1.10) 



Here 



U{Xi,X2 



We solve algebraic system ( |l.9D by using the Newton method according to which cor- 
rection of an approximate solution at each iteration step is performed by the scheme 



g X 



(A: = 1,2,...), 



(1.11) 
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where H{x) is the symmetric Jacobi matrix of functions gi, g2, ■ ■ ■ ,52m of variables xi, 
X2, ■ ■ ■ , X2m, whose elements are calculated as follows 




dui duj , „ , dr 1 

V D\2—UiUj >sds, 

(JjO (JjO kXiO I 



+- 



Ai^Ai 
ldrdU{\i,\2) d 



dvi dvj 
ds ds 



s ds 8X2 



ds 



(ViVj) + 



dU{X2,Xi) \ 

— — h U{X2,Xi) + Qs-r 

XooXo ds 



{j = l,m; i> j), 



sds, 



H,m+j 



+ 



r dU{Xi, A2) dz dr dui dv 



3 



XidXi 


dU{Xi,X2)dz 



ds ds ds ds 

du. 



sdXo 



ds 



+ QX 



dz 



--—Vj — DX2—UiVn > sds 
ds ds ' 



,m). 



So starting from the potential energy stationarity principle, we reduce the problem of 
determining finite deformations of the circular membrane to calculating integrals forming 
algebraic system and elements of the Jacobi matrix with subsequent solving problems of 
linear algebra ( [L.ll ) at each step of the Newton iteration process. The algorithm may 
be effective if one chooses properly systems of coordinate functions {ukis)} and {vk{s)} 
that would allow one to get solutions with high precision, provided the dimension of the 
nonlinear algebraic system is not large. 

To this end, let us expand solutions to be found in a Taylor series in a neighbourhood 
of the point s = 0, using sequential differentiation of equilibrium equations and geometric 
relationships of the shell in order to get coefficients of the expansion. With regard to 
symmetry conditions at the pole of the deformed shell 

Ai = A2 = A, Ti = T2 = T, ki = k2 = k (s = 0), 

we can show that solutions boundered under s — > have the following structure: 



z{s) = ai + a2S^ + 035^ + • ' 
r{s) = bis + b2S^ + bss^ + ■ 



(1.12) 



Taking into account expansions ( 1.12| ), sequences of coordinate functions {uk{s)} and 
{vkis)} satisfying boundary conditions (|1.8D take the form 



Ukis) = is' - l)s'' 



2k~2 



Vk{s) = is' - l)s'' 



2k -1 



{k = l,...,m). 



(1.13) 



Let us present some results of calculations according to the algorithm proposed. To 
illustrate convergence of the Rietz method, we write out values of functions z{s) and 
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r(s) and their first two derivatives, depending on the number m of approximations in 
expansions ( |1.7D at the point s = 0.2 for C = 1.7, D = in Table 1. Ratios of elastic 
constants in Biderman's potential are chosen as follows: 

Ti = 0.02, Ta = -0.015, = 0.00025. 

In the last column of the table we give the relative error of the obtained approximation 
to the solution of the second equilibrium equation 

A = \kiTi + k2T2 - Q\/\C\. (1.14) 



Table 1 



m 


z 


r 


-z' 


r' 


-z" 


-r" 


A 


1 


0.7016 


0.2865 


0.2923 


1.3966 


1.4617 


0.5408 


2 


IQ- 


-1 


2 


0.7824 


0.3026 


0.3975 


1.4644 


1.9279 


0.7252 


6 


10- 


-2 


3 


0.7913 


0.3064 


0.4255 


1.4753 


2.0248 


0.8345 


1 


10' 


-2 


4 


0.7926 


0.3068 


0.4350 


1.4759 


2.0421 


0.8560 


1 


10- 


-3 


5 


0.7926 


0.3069 


0.4362 


1.4758 


2.0415 


0.8591 


6 


10" 


-4 


6 


0.7926 


0.3069 


0.4362 


1.4757 


2.0406 


0.8596 


2 


10" 


-5 



Solutions possess the similar convergence at other points of the interval of integration 
for these equations. Lengthening of the generatrix of the deformed shell affects weakly on 
the degree of convergence of the variational method. For a chosen system of coordinate 
functions, solutions have the uniform convergence as well as their first two derivatives for 
sufficiently wide range of parameters. 

Graphic dependence of the sag of the shell at its pole / = z{0), upon the applied 
uniform pressure {D = 0) is shown in Figure 1. The characteristic indicates the possibility 
of existence of several equilibrium forms for the same value of the load. It is clearly 
seen that the diagram "load-sag" for the material under consideration has two extreme 
points A and B. If the parameter of the load increases monotonely and reaches the value 
of the upper critical load, then there will be an abrupt passage of the system into a new 
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equilibrium state, after which sags of the sheh will again increase smoothly. The monotone 
decrease of the load C (unloading stage) leads to the backward distortion of the shell to 
the initial equilibrium form at lower critical load. Equilibrium states described by a falling 
part of the curve from the point A to the point B are unstable and they are not realized 
at all in the way of loading under consideration. 

In conclusion, we consider to the main principles of choosing an initial approximation of 
the Newton iteration process that allows us to get solutions of nonlinear algebraic system 
(|1.9| ). If solutions of these equations are far from singular points A and B, then one can 
construct an initial approximation, taking advantage of rapid convergence of solutions. As 
we can see from Table 1, the first Rietz approximation is a good solution to the initial 
problem. In this case, we have to solve a system of two nonlinear algebraic equations. 
The initial approximation is determined from physically noncontradictory values of rri 
and X2 (to the positive load there corresponds the positive sag of the shell). Then we use 
the obtained solution to construct P-th approximation, setting the first and (P + l)-th 
components of the vector x equal to obtained values of xi and X2, and putting other 
components equal to zero. The initial approximation constructed in this way improved by 
employing the Newton method with 4-5 iterations. In the case, when it is necessary to 
obtain a series of equilibrium states, when choosing initial values of the vector x one can 
use different extrapolational formulae to prolong solutions in the parameter C. However, 
when implementing the procedure of prolonging a solution with approximation to critical 
points A and B, the convergence of the Newton process sharply comes down, and then 
the divergence appears clearly. To go through these equilibrium states of the shell one 
should change prolongation parameters. The main point of this approach is that starting 
from some moment, the parameter of the load is supposed to be unknown, and as a given 
parameter we choose the value of the sag of the shell at the pole. Artificial introduction of 
a new parameter provides that the uniqueness condition for a solution in a neighbourhood 
of a fixed parameter is fulfilled, and, therefore, Jacobian is nondegenerate, which is a 
necessary condition for the Newton method to be applicable. 

The algorithm proposed is also effective for other rotation shells in the nondeformed 
state, the generatrix of which intersects the symmetry axis at the right angle and principle 
curvatures at its pole are equal. 

2 Adaptive Rietz method under conditions of singular per- 
turbance of the initial problem 

Let us consider the case of loading a membrane by a hydrostatic load for C* <C 1 and 
D ^ 1 (heavy liquid). For such parameters of the load, stresses in the median surface of 
the deformed shell will be small, i.e., {Ti,T2} <C 1, and the shell will take a nearly plane 
form with the exception of a narrow domain near the fixed contour of the membrane. 
The increase of the constant component of the load C* leads to the increase of stresses 
in the membrane, that begin to play the prevalent role in forming its deformed surface. 
This is accompanied by increasing the width of a domain with rapidly changing functions 
of movements and stresses, and turning a plane part of the surface of the shell into a 
curvilinear surface. 

After introducing a small parameter /x = 1/D, equilibrium equations (|1.6| ) in the coor- 



Variational Methods for Solving Nonlinear Boundary Problems of Statics 



41 



dinate system associated with the center of the nondeformed membrane can be rewritten 
as follows: 



Li{s,r,z',r',z",r")=0, 
fi[L2{s,r,z',r',z",r")-C*]=l-z, 



(2.1) 



where 



dT 1 dv 

Li = -^ + --^(Ti - Ts), La = hTi + fcsTs, / = z(0). 
as r as 

In the case under consideration, the presence of a small parameter at higher derivatives 
transfers the initial regular problem into a class of nonlinear singular perturbed problems. 
In this case, solutions will have domains both with fast and slow variability. As a result, 
when approximating solutions by means of polynomial basis, the convergence of successive 
approximations in the Rietz method sharply comes down, and the increase of the number 



of coordinate functions in expansions (O) leads to a loss of calculation stability before 
the required precision is achieved. The situation is also aggravate by that the singularity 
degree of the problem depends on a solution itself that is a priori unknown. Thus, there 
arise the principal problem of extending a class of admissible functions by functions that 
allow getting within the framework of a unique algorithm approximate solutions with the 
equally high precision both in the presence of a narrow area of great changes in solutions 
and without it. 

The problem indicated plays a key role in realizing the method of variations when 
solving singularly perturbed boundary problems. It is one of suffuciently difficult problems 
of mathematical physics and theory of function approximation. Apparently for this reason, 
variantional methods (unlike finite difference methods |^]) are not developed properly for 
solving even linear boundary problems with a small parameter the higher derivative. 

One of the most important merits of the method of variations is that when constructing 
approximate solutions, one can use an information about problem, that can be obtained 
by means of preliminary analysis of required solutions. On the other hand, methods of a 
small parameter may be also used to determine a structure of a solution and characterize 
its degeneration with the parameter at the higher derivative tending to zero. With this 
objective in view, let us carry out an asymptotic analysis of the initial equations. 

Searching for a direct expansion of solutions to system (|2.lD in the form 



z{s) = X^/zV'=)(s), r{s) = ^ A^'H^), (2.2) 

k=0 k=0 

let us expand operators Lj (i = 1, 2) in a Taylor series in the parameter /i 



Having substituted (2^) and ( p.3[ ) into ( |2.1| ) and having equated coefficients of expansions 
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of the same degrees of /x, we come to 

'^^'^ 0, (L2),=o-C* = -z«(.), 



^/^ / M=o (2.4) 
"'^ 0. ^fS) ^ -.'")(.) („^2.3....). 

M=0 



d/i" y ' (n-1)! I 



We can obtain a bounded solution of the nonhnear equation {J-i\)p,=Q = with respect to 
the function A^'^ (s) by expanding it in a Taylor series in a neighbourhood of the point ,3 = 
with regard to the symmetry at a pole of the deformed shell. In this case, we can show 
that 

r(°)(s) = ais, z(o)(s) = Z, (2.5) 

where ai is an arbitrary constant. 

Furthermore, with regard for expressions ( |2.5D , linear equations of the first approxima- 
tion allow us to get expressions for functions r^)(s) and z^^\s): 

r«(s) = a2S, zW(s) = C7*. 

Having performed several transformations that are simple but sufficiently awkward, 
we can show that 

A^\s) = anS, z^^Hs)=0 (n = 2,3,...), 

where a„ are arbitrary constants. 

Collecting the results obtained, we come to the general form of a solution bounded for 
s = as a result of direct expansion (p. 21): 



r{s) = A{fi)s, z{s) = C/D, C = C* + ID, (2.6) 

where A{^) is a function of /x; C is the constant component of the hydrostatic pressure, 
that differs from C* at the expense of that the origin of the coordinate system passes into 
the center of the nondeformed membrane. 



Proceeding from equations (2.1), it is not possible to establish the structure of solutions 
in the boundary area, since equations of the first approximation are nonlinear and they 
are not integrated in the explicit form. For this reason, in what follows, using the results 
of the asymptotic analysis of the linearized equilibrium equations with respect to some 



equilibrium state that is described by an analytic function of the form (1.12), we propose to 
establish the structure of their solutions and after that to endow solutions of the nonlinear 
problem with these properties. 

Suppose that the membrane is in the deformed state as a result of applying the hy- 
drostatic load Q = C* — Dz. Let the shell have passed into a new state at the expense 
of changing the constant pressure component by a small value 5C* . Denote by u and w 
projections of a movement of the shell from this load onto directions of its generatrix and 
outer normal, respectively. We shall characterize the main strained deformed state by 
stresses Ti and T2, and principal degrees of lengthenings Ai, A2 and A3. Having performed 
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linearization of elasticity relations ( |l.lD , we can establish the following relation between 
additional stresses STi and 6T2 and deformations of the shell: 

5Ti = ciiEi + 01262, ST2 = C2iei + C22e2, (2.7) 

where 

1 dti 1 dr 

ei = —— + kiw, 82 = ^--^u + k2w; 

Ai as r\i as 

cii = /i(Ai, A2), C12 = /2(Ai, A2), C21 = /2(A2, Ai), C22 = /i(A2, Ai); 



/i(Ai, A2) = (A? A3 + 3Ai) ^ + (A1A2 + 3XlXf 



+2A3 (A? - Xl) {Au + 2^12A2 + ^22At) , 
/2(Ai, A2) = (3Ai - A? As) ^ + (A1A2 + X^Xl) ^ 

+2A3 (A? - Xl) {Xl - Ai) (^11 + Au {Xl + Xl) + A22XIXI) , 

^ifc = 7T777r {i,k = 1,2). 
olid Ik 

Linearized equilibrium equations in movements, describing the perturbed state of the 
shell are of the form 



Lii(n) + L12H = 0, 

—w + —u \ = {e^dC* + AI) rXi 



(2.8) 



where 



T / \ d f du\ T / \ dw 

Lii(u) = -— \ai— ] + a2U, Luiw) = 03— + 04^, 
as \ as J as 

L22{w) = --7- ( Pi^] + P2W, L2i{u) = + P4U, 

ds \ ds J ds 

Cur (1) daf^ (i) C22 ( dr\'^ 



(2) 
4 



(2) C12 + Ti (ir / ; I / ^ (1) .^a 
02 = ^ 3-, as = -?"(C11K1 + C2lfc2), 04 = 0^4 H 1— 

Ai ds ds 
a^i' = [(c2i + T2)ki + 022/^2]^, ^ = -?^[ciiA;i + (ci2 + ^1)^2], 

/?i = /?2 = -rAi [(cii - Ti)/;;! + (ci2 + C2i)A;i/c2 + (c22 - 72)fei] ; 
Ai 

= -as, Pa = rf^ + rf^ = [cisfei + (C22 - ^2)^:2]^, 

/3f)=rTiA:i, = /U, A/ = u'(O). 
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One can show that the structure of a direct expansion of solutions to equations in 
the parameter e corresponds to structure ( |2.6| ) estabUshed for the system of initial nonlin- 
ear equations. The subject of further research is to construct integrals with high variability 
of a system of linear equations with variable coefficients ( |2.8D . Since a direct expansion is 
an approximation to some solution of initial equations ( p.8D , we shall seek boundary layer 
part of an asymptotics as an approximation to a solution of the relative homogeneous 
system of equations. In so doing, we shall take into account that functions r{s) and z(s) 
characterizing the geometry of the shell in the main deformed state are analytic functions 
in a certain domain and according to ( 1.12| ) they are expanded in power series in even 



and odd powers of s, respectively. On the basis of this fact, we can rewrite the initial 
homogeneous system of equations as follows: 



2 d?u du dw 

s -j-iT + spn— + spi2— + qnu + Quw = 0, 
as^ ds ds 

2 od W 2 1 "'"^ \ 

{sP21-J^ + SP22-J^j+q2lU + q22W = 0, 

where Pij{s) and qij{s) are analytic functions for s € [0, 1]. 

Proceeding from the form of system of equations (|2.9D , we conclude that this system 
has singularities both in the parameter and unknown variable. The presence of a small 
parameter at the higher derivative results in existence of a narrow area in a neighbourhood 
of the point s = 1, in which solutions of the system have large gradients. The point s = 0, 
in its turn, is a regular point for equations ( |2.9D , that assumes the relative asymptotics 
of solutions when s — s- 0. Since our final goal is to construct an algorithm for solving a 
problem with the uniform convergence in the parameter e, in what follows it is necessary 
to take into account both of these facts. 

We have not succeededin constructing solutions a singular perturbed equation, with a 
regular singular point, that belongs to the class of equations under consideration containing 
an exponent factor [pX|]. Simplified equations obtained from the initial ones while only the 
first terms of expansion of their coefficients remain unchanged, may serve as a guiding 
line for an asymptotic representation of solutions to system ( |2.9| ). In this case, we can 
establish that a regular solution of these equations is expressed by Bessel functions of 
order zero and their first derivatives. According to what has been said, we seek solutions 
of homogeneous system of equations ( |2.8| ) with a large variability index in the following 
form: 

u = ui{s,e)Jo{x) + U2(s,e)Jg(j;), w = wi{s,e)Jo{x) + W2{s , e) Jq{x) , 

T f (2-10) 

X = -, r = / P{t)dt. 



Here Ui{s,e), 'Wi{s,e) {i = 1,2) and P{t) are functions to be determined. 

The advisability of representing solutions in the form ( |2.1C1| ) will be warranted if func- 
tions Ui{s,£) and t(;j(s,e) that we obtain are regular and can be presented as a direct 
expansion in the parameter e. 

Having substituted expressions ( |2.10| ) into homogeneous system of equations ( |2.8D and 
set equal to zero coefficients of Joix) and Jq{x), we get four equations for the functions 
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Ui{s,e), Wi{s,e) and P{s): 

Lii(ni) + L22{wi) + {aiPu2y/e + (uiP^a^/e^ 

+ai (^u'^P - /e - {a^W2P)/e = 0, 

-^^11(^2) + Li2{w2) + {aiPu2/T)' - (aiPuiY/e + {a3Pwi)/e - {a3Pw2)/T 

+aiU2P'^{l/e^ - l/r^) + Pai(niP/r - u[)/e + n'aPai/r = 0, ^2.11) 
[L2i(ti2) + L22{w2) + (/?iW2 P/t)' + P{PlW2 - P3U2)/t] + rW2r' + rU2Z 
+e [(53U1P - iPiWiPy + (3iP{wiP/T - w[)] + f3iP^W2{l - eV^') = 0, 

e^[L2i{ui) + L22K)] + e [{l3iW2P)' + /3iP(u;^ - W2P/t) - 133U2P] 
+wiP'^Pi + rwir' + ruiz' = 0. 



We shall seek a solution to system of equations ( 2.11 ) in the form 



00 00 

ui — 2_^e {s), U2 — 2_^e U2 {s) 

k=0 k=0 

00 00 



(2.12) 



k=0 k=0 

To determine functions uf'\ w^^^ and P{s) we substitute expansions (|2.12|) into equa- 
tions ( |2.11| ) and set coefficients of the same degrees of e equal to zero. 

In the zeroth-order approximation, i.e., for functions n^'^'* and w^^\ we have 

P^aiuf ^ - asPw;^^^ = V, Pai4°^ + a^wf^ = 0, (p2/3i + rr')wf'^ = 0, 



(0) ^ _a34^ „2 ^ ^ _ Afcosa 

2 aiP ' I3ids Ti ' ^ ' ' 

Here a is the angle formed by the outer normal to the deformed shell and its symmetry 
axis. 



It follows from (|2l^ ) that P{s) is an even and pure imaginary function for < a < 



The fourth equation with regard for ( ^.13| ) takes the form 

^dw\^^ (P 1 „N / ^^ (0) / N a^dz Idr 
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Hence up to an arbitrary constant, we obtain a regular solution for function wf^^: 



= y ^^^P J 'Pi^)ds- (2.15) 

Therefore, after the zeroth-order approximation, the function remains a free func- 
tion that win be determined at the second step of the asymptotic integration of equa- 
tions (|1|). 



Without to finding exphcit expressions for terms in expansions (2.12) of higher ap- 



proximations, one may conclude from the form of system ( |2.11| ) with regard to expressions 
for its coefficients, that u^\s) and w[^\s) are even functions, and uf'\s) and w^\s) 
are odd ones. Specific expressions for these functions after the zeroth-order and first 
approximations show that they are regular for s = 0. 

Putting together the results obtained, we conclude that a solution to system ( |2.8D , 
bounded for s = and having the large variability, can be presented in the form of the 
following formal expansions: 

oo oo 



k=0 k=0 

2.16 

oo oo 

wis, e) = h{y) X: e^'wf^ (s) + /^(y) J2 '"'"''^^'^ 



k=0 k=0 



where Iq is a Bessel function of the order zero of pure imaginary argument; y is a function 
regular and odd in argument s. 

The constructed asymptotic representations of solutions (2^) and ( 2.16| ) quantitatively 



present the structure of a solution, that we use in a sequel to construct systems of coor- 
dinate functions when solving the initial problem by the method of variations. 

Taking into account the asymptotics of expansions ( |2.16| ) for s ^ 0, and also the 
asymptotics of fucntions Inix) for large values of the argument, after passing from normal 
and tangent movements of the shell to its movements in axial and radial directions, the 
general form of solutions to the initial nonlinear problem can be presented as follows: 



z{s) = ao + ipp{s){ai + a2S^ + a^s'^ H ), 

r{s) = bos + ipp{s){bis + 62^^ + ^3^^ -\ ), 

\fc=0 / 



(2.17) 



where Oj, hi and pi are arbitrary constants. 

Subjecting expressions ( 2.17| ) to boundary conditions (1^), coordinate systems {uk{s)} 
and {vk{s)} take the form 



Ul{s,Pi) = [l — , U2{s,Pi) = — -— (s -1), 

V'p(l)/ '/'p(l) (2.18) 



Uk{s,pi) = s^Uk^i{s,pi) (A; = 3,4, . . .), Vk{s,pi) = suk{s,pi). 
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If constants pi are known, then solving nonlinear algebraic system (|l.9|) by means of 
the iteration Newton method, we can construct an approximate solution for the problem. 
However, the question of specific values of parameters pi for the nonlinear boundary 
problem under consideration remains open. Hence, from the conditions dl/dpi = 0, 
z = l,2,...,nwe can get n additional equations with respect to parameters pi. These 
equations are of the following form: 



' ds ds \ dpi J ds dpi ' ds ds \ dpi 



(2.19) 



+ (f/(A„M- + QA.-j - 



sds. 



It is convenient to look for a solution of system (2.19) by the generalized secant 
method [|l^]. The whole algorithm is reduced to that at each step of changing param- 
eters Pi (outer iteration process) one needs to solve nonlinear algebraic system (|1.9D. 

Therefore, unlike the traditional Rietz method of solving nonlinear boundary problems, 
a peculiarity of the given approach is that at the expense of choosing parameters pi that 
characterize the variability of required solutions, by using a computer, we construct a 
system of coordinate functions, that approximates solutions to the initial problem in a 
way optimal in some sense. 

In conclusion, let us present results of calculations, giving an idea of possibilities and 
effectiveness of the algorithm proposed. To illustrate the convergence of successive ap- 
proximations, in Table 2 we write out values of functions z{s) and r{s) and their first 
two derivatives depending on the number m of terms in expansions (|l]^) at the point 
s = 0.9. In the last column we give the relative error of an approximate solution for 
problem ( 1.14| ). Parameters of the hydrostatic load in the coordinate system associated 
with the center of the nondeformed membrane are C = 0.5, D = 10. Ratios of constants 
in elastic potential (|1.2|) are chosen as follows: 



0.1, 



r. 



0. 



All calculations are performed with regard for only one parameter Pi in coordinate systems 

(HI). 



Table 2 



m 


z- 10 


-z' 


-z" 


r 


r' 


-r" 


A 


1 


0.32896 


0.18816 


1.9342 


0.90532 


0.97550 


0.38008 


4 


10" 


-1 


2 


0.36718 


0.18415 


2.4584 


0.90694 


0.99159 


0.46372 


1 


10- 


-2 


3 


0.36533 


0.18145 


2.3754 


0.90692 


0.99250 


0.42182 


1 


10" 


-2 


4 


0.36451 


0.17825 


2.3497 


0.90693 


0.99278 


0.41433 


3 


10" 


-4 


5 


0.36450 


0.17830 


2.3498 


0.90693 


0.99276 


0.41456 


2 


10- 


-4 


6 


0.36448 


0.17841 


2.3461 


0.90693 


0.99275 


0.41404 


3 


10" 


-5 



Solutions possess also the analogous convergence at other points from the interval of 
integration of the initial equations. Here the high precision of an approximate solution is 
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achieved at the expense of both increasing the dimension of algebraic system ( |l.9D and 
choosing parameters of the optimal coordinate system. It turns out that it is sufficient 
to restrict oneself to one parameter pi only. Calculations show that requirements to the 
precision of determining parameters pi concern with input data of the hydrostatic load 
with high singularity in solutions of the problem only. As the width of the boundary 
layer increases, the rate of convergence of the algorithm is preserved, but influence of the 
precision of determining parameters for the optimal coordinate system on the final result 
decreases. There comes a moment when for some value of the hydrostatic load we do not 
need any more to organize the outer iteration process for seeking for values of parameters 
Pi. In this case, the rate of convergence of the given algorithm coincides with the rate of 
convergence of the algorithm based on the use of an exponent basis. It is worth noting 
that a polynomial basis for the example considered above allows one to obtain the relative 
error A for the first eight terms in expansions ( |1.7D , that is equal to ~ 10~^. The further 
increase of the dimension of algebraic system ( |l.9D in this case leads to the loss of the 
stability of calculations. 



While performing the calculating experiment, it came out that functional (L4) in the 
space of parameters pi had only one minimum. This means that for the approximation 
under consideration there exists only one optimal coordinate system, which considerably 
simplify its search by means of a computer. 

Therefore, the proposed approach allows one to obtain approximate solutions with a 
suffuciently high degree of precision for the small dimension of an algebraic system both 
in the presence of the singular perturbation of a problem and without it. 

Meridional cuts of the deformed shell for different values of the hydrostatic load C 
are presented in Figure 2. The other input data coincide with those used for calculations 
presented in Table 2. Dashed lines correspond to paths of passage of points of the non- 
deformed membrane to points of its deformed surface. Curves of dependence of principal 
degrees of lengthenings Ai and stresses Ti (continuous curves) and also A2 and stresses T2 
(dashed curves) upon the parameter s for different values of the parameter C are presented 
in Figures 3 and 4. From figures one can see that at the starting stage of deforming the 
membrane, the main part of its surface is nearly plane that is in the homogeneous biaxial 
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Figure 3 Figure 4 

strained state. Only in a neighbourhood of the support contour of the shell its surface 
and intererior stresses undergo considerable changes. As the constant component of the 
hydrostatic load C increases, the plate part of the surface gradually decreases. 

Since it is convenient to investigate behaviour of material of elastomers for large de- 
formations by means of experiments on homogeneous biaxial extension of a thin sheet 
(that often causes some difficulties when realizing such a type of deformations for ex- 
perimental tests), the presence of the homogeneous strained state in the central part of 
the deformed membrane under the load of heavy liquid may be used in establishing a 
functional dependence of an elastic potential upon deformation invariants. 
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